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We present the first measurements at a hadron collider of differential cross sections for 
Z/y*+jet+X production in A<p(Z,]et), |Aj/(Z,jet)| and \yboost(Z + jet) |. Vector boson produc- 
tion in association with jets is an excellent probe of QCD and constitutes the main background to 
many small cross section processes, such as associated Higgs production. These measurements are 
crucial tests of the predictions of perturbative QCD and current event generators, which have varied 
success in describing the data. Using these measurements as inputs in tuning event generators will 
increase the experimental sensitivity to rare signals. 

PACS numbers: 12.38.Qk, 13.85.Qk, 13.87.-a 



The production of the massive vector bosons W and 
Z at hadron colliders, such as the Fermilab Tevatron 
and CERN Large Hadron Collider, has distinctive fea- 
tures which lead to many applications. The electron and 
muon decay modes of these bosons are experimentally 



rather simple to identify in the complex environment 
of a hadron collider and typically have very little back- 
ground. Reconstruction of the boson kinematics from 
its decay products provides a unique, colorless probe of 
the underlying hadron collision and any hadronic recoil 
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to the boson. The production of high energy jets in as- 
sociation with W and Z bosons (V+jets) typically lies 
within the regime of perturbative quantum chromody- 
namics (pQCD), however such final states are difficult to 
calculate to higher orders in perturbation theory. Several 
tools have been developed for generating V+jets events, 
from the parton-shower event generators pythia [l[ and 
herwig [2|, to generators that combine tree level ma- 
trix element calculations with parton showers, such as 
ALPGEN 3 (using pythia or HERWIG for showering and 
hadronization) and SHERPA Comparisons between 
these generators show they suffer from significant uncer- 
tainties and differ in the predicted kinematics of V+jet 
production Studies of such final states are therefore 
an excellent testing ground for theoretical predictions, 
and inputs from measurements are needed to improve 
these models. V+jets is also an experimental signature 
of many other processes with significantly smaller cross 
sections, such as top quark pair production, associated 
production of the Higgs boson with a W or Z boson, 
and the decays of particles within many supersymmetric 
scenarios. Identifying V+jets final states resulting from 
such rare processes relies upon a precise understanding 
of the more copious V+jet production predicted purely 
by QCD. 

Previous measurements at the Tevatron have studied 
the kinematics of inclusive Zj 7* production [1, 0, 0] , of 
the Z/j* in events with at least one jet M, of the jets 
in Z/7* and W events d EE EI El El] , and of the 
production of ZH* and W in association with heavy fla- 



vor 



QUEUE 

I17H . In this Letter, we describe the first 
measurements of the angular correlations between the 
Z/7* and leading jet in Z/7*+jct+A" production. Using 
the decay mode Z/j* — > /i/i, differential Z/j*+jet+X 
cross sections are measured, binned in the azimuthal an- 
gle between the Z/j* and leading jet, A^Z, jet), the 
absolute value of the rapidity [18| difference between 
the Z/j* and leading jet, \Ay(Z, jet) | , and the absolute 
value of the average rapidity of the Z/-f* and leading 
jet, \yboost{Z + jet) I . These differential cross sections are 
normalized to the measured inclusive Z/j* cross section, 
cancelling many systematic uncertainties. As for a previ- 
ous measurement of A0(jet, jet) in inclusive two jet pro- 
duction [ill, the A</>(Z,jet) distribution is sensitive to 
QCD radiation. In the absence of additional radiation, 
the Z/"f* and jet would be produced with equal and op- 
posite transverse momenta, with A</>(Z, jet)= n. Further 
non-collinear radiation at any transverse momentum, px, 
results in A0(Z,jet) deviating from 7r, with higher px 
radiation giving a smaller A</>(Z,jet). The rapidity vari- 
ables, Ay(Z,jet) and ?/boost(Z + jet), have primary con- 
tributions from the relative momenta of the incoming 
partons in the hard scatter, but again are modified by 
any additional QCD radiation. These measurements are 
therefore excellent tests of the inclusion of QCD radiation 
in theoretical models, without requiring that more jets 



be observed in the event. As a result, the measurements 
avoid the experimental uncertainties that would be asso- 
ciated with requiring additional jets to be reconstructed, 
and is also sensitive to jets below detector reconstruction 
thresholds. 

This analysis uses a dataset of pp collisions at 
yfs = 1.96 TeV, corresponding to an integrated luminos- 
ity of 0.97 ± 0.06 ib" 1 [20I recorded by the DO detector 
between April 2002 and February 2006. A full description 
of the DO detector is available elsewhere [2l[ , and only the 
components most relevant to this analysis are described 
here. Immediately surrounding the pp interaction region 
are two tracking detectors: a silicon microstrip tracker 
and a scintillating fiber tracker, housed inside a solcnoidal 
magnet providing a field of approximately 2 T. These 
trackers are used to measure the momenta of charged 
particles and to reconstruct the primary interaction point 
in each collision. Outside the solenoid lies a liquid argon 
and uranium calorimeter which is split into three sec- 
tions: a central section extending to \rj\ < 1.1 [22| and 
two forward sections covering 1.4 < \rj\ < 4. Scintillat- 
ing detectors provides additional energy measurements 
between the central and forward calorimeters. Outside 
the calorimeter lie three layers of muon detectors which 
are a combination of scintillating pixels and drift tubes. 
Between the first and second layers lies a 1.8 T iron 
toroidal magnet, providing an independent muon mo- 
mentum measurement used in the trigger system. 

Events used in this analysis are selected by at least 
one of a suite of single-muon triggers. Each of these trig- 
gers uses fast readout from the muon system scintilla- 
tors and fiber tracker to initially identify events, then 
information from the full tracking and muon systems 
to provide further rejection. Additional requirements 
are then applied to obtain a sample of Z/j* candidate 
events. Using information from the muon detectors and 
the tracking system, two muons of opposite charge and 
Pt> 15 GeV are required, with a dimuon mass in the 
range 65 < < 115 GeV. To reject cosmic rays 

and poorly reconstructed muons, the muon tracks are 
required to match the reconstructed primary interaction 
point both transverse and parallel to the beam direction; 
the two muon tracks are also required not to be collinear. 
Finally, the muons are required to be consistent with the 
pp bunch crossing time, using timing information from 
the muon system scintillators. 

Jets are reconstructed using the DO Run II seeded, iter- 
ative mid-point cone algorithm [23| on clusters of energy 
deposited in the calorimeter. The algorithm is config- 
ured with a split-merge fraction of 0.5 and cone radius 
of y/ (A(f>) 2 + (Ay) 2 — 0.5. Shape and quality cuts re- 
ject spurious jets caused by electrons, photons or noise 
in the calorimeter. Further corrections are applied for 
the calorimeter response, instrumental out-of-cone show- 
ering effects, and additional energy deposits caused by 
instrumental noise and pile-up from multiple pp interac- 
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tions and previous pp bunch crossings. These corrections 
are derived by balancing the px in 7 + jet events, where 
the 7 and jet are opposite in </>. After corrections, jets 
with p T > 20 GeV are selected. 

Further selections are applied to limit the measure- 
ment to regions with high detection efficiency and well- 
understood detector performance: the muons are re- 
quired to have |ry| < 1.7, the primary vertex must lie 
within 50 cm of the center of the detector along the di- 
rection of the beam, and only jets with \y\ < 2.8 are 
considered. A total of 59,336 Z/j* — ► candidate 
events are selected before jet requirements, of which 9,927 
contain at least one jet passing all selections. In events 
containing more than one jet, the highest px jet is se- 
lected to calculate the angular variables. Finally, the low 
p^ region is excluded as at low p^, the measurement of 
the Z/j* azimuthal angle is dominated by the experi- 
mental resolution of the muon px- Additionally, these 
events may contain a significant fraction of jets from ad- 
ditional parton interactions, essentially uncorrelated with 
the Z/j* production. For these reasons, the measure- 
ment is carried out in two kinematic regions: first, with 
Px> 25 GeV, and then also raising the pf selection to 
> 45 GeV to probe the higher px region in more detail. 
These contain 5900 and 2449 events respectively. 

The main source of background in this analysis is 
muons from semi-lcptonic decays in high energy jets or 
VF+jet production. This is reduced to negligible levels 
by limiting the sum of track px and the calorimeter en- 
ergy allowed in a cone around each muon, and avoiding 
the overlap between muons and jets by requiring angular 
separation AR = y/ (A(f>) 2 + (Arj) 2 > 0.5. The remain- 
ing contribution of M^+jets is extracted from data by 
studying the distribution of muons failing these require- 
ments and extrapolating into the signal region, and is 
found to be < 0.5% of the final sample. Other sources 
of background (top quark, diboson, and Z/j* — > t + t~ 
production) are estimated using pythia simulation us- 
ing CTEQ6L1 parton distribution functions (PDFs) [H, 
and the "Tune- A" underlying event settings [jjij ]. passed 
through a GEANT-based [26[ simulation of the DO detec- 
tor, and normalized to higher order theoretical predic- 



tions 
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28J. A 10% systematic uncertainty is assigned 
to this normalization. Backgrounds are found to be neg- 
ligible almost everywhere, however the top quark contri- 
bution is evenly distributed in A<p(Z, jet) and accounts 
for up to 11% of the data at small values of A<j)(Z, jet). 
The estimated background contributions are subtracted 
from all data bins. 

To extract differential cross sections, the measured 
events must be corrected to the particle level [2j|, ac- 
counting for detector resolution, acceptance, and effi- 
ciency. These corrections are derived from simulated 
Z/7*+jet events generated with ALPGEN v2.11 using 
CTEQ6L1 PDFs, and showered using pythia v6.413 
with the "Tune- A" underlying event settings. These sim- 



ulated events are then passed through a GEANT-based 
simulation of the DO detector, and real data events from 
random bunch crossings are overlaid to reproduce the ef- 
fects of multiple pp interactions and detector noise. 

Further corrections are applied to the simulation in 
order to improve the description of the data. The muon 
trigger is not simulated; instead the trigger efficiency is 
measured in data, and parameterized in terms of the ge- 
ometry of the muon system. This efficiency is then ap- 
plied on an event-by-event basis to the simulation, with 
the average efficiency being approximately 88%. The ef- 
ficiencies of the muon reconstruction and isolation re- 
quirements are measured in data and in the simulations, 
and adjustments of order 3% are applied to the simula- 
tion to correct for the differences. The muon px resolu- 
tion is studied by comparing the shape of the Z/j* mass 
peak in data and the simulation, and further smearing on 
1/pt of order 5% is applied to the simulation to repro- 
duce the data. Jet corrections are studied by measuring 
the px balance in back-to-back Z/7*+jet configurations, 
and further scaling and smearing, which reaches a max- 
imum of 4% at low px, is applied to the simulation to 
match data. Finally, kinematic variables as produced by 
the event generator are re-weighted to provide a good 
simultaneous description of the variables important to 
this analysis: pf , y z , p J ^*, y jot , as well as A</>(Z,jet), 
Ay(Z, jet)and yboost(-Z' + jet). After these corrections, 
the simulation provides a good description of all mea- 
sured distributions. 

To then correct for detector effects, the corresponding 
particle level quantities must first be defined. To min- 
imize dependence on models of particle production, de- 
cay and radiation, the definitions are based on the stable 
particles that enter the detector in the simulation (with- 
out identifying the origin of these particles) and corre- 
spond as closely as possible to the detector level defini- 
tions of muons and jets. The Z/j* candidate is recon- 
structed from the highest-p^ P + and fi~ with |y M | < 1.7, 
then requiring 65 < M MM < 115 GeV. As for the detec- 
tor level analysis, two samples with pf > 25 GeV and 
Px> 45 GeV are selected. The two muons used to recon- 
struct the Z/'y* and any photons in a cone of AR < 0.2 
around those muons (dominated by QED radiation from 
the muons), are excluded from the jet reconstruction. All 
other stable particles are passed to the DO Run II jet al- 
gorithm, with the same settings as for detector level jets. 
Jets with \y\ < 2.8 are considered, then the highest px 
jet is selected and required to have px > 20 GeV. 

As the distributions of Ay(ZJet) and ?/boost(Z + jet) 
are symmetric around zero, the absolute value of each 
is considered in order to increase the statistical preci- 
sion. The distributions of A0(Z,jet), |Ay(Z,jet)|, and 

|2/boost 

(Z+jet)| are then binned, with the binning de- 
termined by a combination of detector resolution consid- 
erations, maintaining reasonable numbers of data events 
in each bin, and maximizing sensitivity to shape differ- 
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ences predicted by different models of Z/j*+}et produc- 
tion. With the chosen binning, the detector resolution 
causes little migration between bins (bin purities are gen- 
erally above 90%, with the purity defined as the fraction 
of events in a given particle level bin which are in the 
same bin at detector level). Purities for A<j>(Z,jet) drop 
to around 70 % in some bins, due to the worse exper- 
imental resolution. However, in all cases the main ef- 
fect to correct is detector acceptance and efficiency. To 
do this, the distributions of A</>(Z,jet), |Ay(Z,jet)|, and 
|?/boost {Z + jet) | are populated in the simulation indepen- 
dently for the particle level and detector level selections. 
The ratio of particle level to detector level is then ap- 
plied to data, with a typical correction factor being ap- 
proximately 2.2, with some dependence on the variable 
under consideration. The statistical uncertainty assigned 
to the resulting differential cross section corresponds to 
the detector level statistical uncertainty. The data are 
then normalized to the total Z/j* production cross sec- 
tion (with no jet or pj, requirements, but the same muon 
\y\ and dimuon mass requirements) measured with this 
sample. 

Finally, systematic uncertainties are assessed. Several 
sources are considered, beginning with the uncertainties 
on pt resolution, energy scale and detector efficiencies 
for both muons and jets. These corrections are shifted 
individually up and down one standard deviation, and 
the ratio of particle level to detector level re-derived. 
The difference in the final result is assigned as a sys- 
tematic uncertainty. The effects of the various kinematic 
re-weightings are assessed by turning off each of these 
in turn and repeating the analysis; however, due to the 
small effects of detector resolution, the results are largely 
insensitive to these tests except for A<f)(Z, jet), where the 
lower bin purities lead to a larger correlation between the 
shape of the distribution and the detector to particle-level 
correction factor. The region of low A0(Z,jet) is found 
to be particularly sensitive to jets from additional inter- 
actions in an event, so a further uncertainty is assessed 
by varying the real data from random bunch crossings 
overlaid on the simulation. The systematic uncertainties 
are combined in quadrature, with the jet energy scale 
uncertainty generally being the largest single source for 
\Ay(Z, jet) | and |yboost(-^ + jet) j , and accounting for ap- 
proximately half of the total uncertainty. For A0(Z,jet), 
the uncertainty on the kinematic re- weightings is compa- 
rable to or larger than the jet energy scale uncertainty. 
For the selection with > 25 GeV, the total systematic 
uncertainty is of comparable size to the statistical uncer- 
tainty; for the selection withp^> 45 GeV, the statistical 
uncertainty dominates. 

Predictions for A<f>(Z, jet), |Ay(Z,jet)| and 
|yt>oost(-^ + jet) | distributions are obtained from several 
theoretical models, as well as leading order (LO) and 
next-to-leading order (NLO) pQCD calculations. Predic- 
tions at LO and NLO are obtained with MCFM v5.6 [27j |. 



together with the MSTW2008 LO and NLO PDFs 
respectively. The distribution at A<fi(Z , jet)=n contains 
divergences and is excluded from comparisons to data. 
Re-normalization and factorization scales are set to the 
sum in quadrature of the mass and pt of the Z/j* in 
each event, and the dependence on this choice is assessed 
by varying both scales simultaneously up and down by 
a factor of two, both for the differential distribution and 
the inclusive Z/j* cross section used in normalization. 
PDF uncertainties are assessed using the MSTW2008 
68% error sets, again taking into account the effect on 
the differential distribution and the inclusive Z/j* cross 
section used in normalization. These are found to be 
approximately a factor of two smaller than the scale 
uncertainties. Comparisons with data are performed 
after correcting the parton level prediction from MCFM 
for the effects of hadronization and the underlying 
event. These corrections have been derived from a 
sample of Z/7*+jet events generated with pythia 
v6.421 [lH using the underlying event tunc QW [HI 
with the CTEQ6.1M PDFs 0. They are derived by 
comparing the full prediction (taken from the final 
state particles, including the underlying event) to the 
purely perturbative part (calculated from partons taken 
after the parton shower, with no underlying event), and 
are typically around 4%. However, the low A(f>(Z,]et) 
(< 1.5 rad) region is dominated by non-perturbative 
effects, so the pQCD calculation for this bin is excluded. 
Corrections for quantum clcctrodynamic final state 
radiation (FSR) from the muons are also derived from 
the same pythia sample, by comparing the prediction 
calculated using the muons after FSR to those using the 
generated boson. These are typically less than 1% after 
accounting for the effect on the inclusive Z/-f* cross 
section and are primarily the result of events migrating 
out of the mass window. 

Predictions are also obtained from four current event 
generators. When considering the number of generators 
available and the various tunes of those generators, a 
complete survey of the field would be impossible. We 
choose to focus on the current matrix element calcula- 
tions with matched parton showers, as implemented in 
SHERPA and ALPGEN, as previous measurements indicate 
these provide the best description of boson+jets final 
states 0, The ALPGEN matrix element calculation 
can be interfaced to the pythia or herwig parton shower 
and hadronization models, and we test both. Further, 
to assess the impact of the additional matrix elements 
in ALPGEN, we also run pythia and herwig stand- 
alone. Unless stated otherwise, we use the same PDF 
set throughout: CTEQ6.1M. First, a sample of events is 
generated with SHERPA vl.1.3 with up to three partons 
in the matrix element calculation. Parton jets from the 
matrix element are required to have px> 13 GeV and 
Ai?(jet, jet) > 0.4. In SHERPA, both the renormalization 
and factorization scales are set according to the CKKW 
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prescription 34|. To provide a typical uncertainty on a 
prediction from an event generator, SHERPA samples are 
also generated with the renormalization and factorization 
scales varied up and down by a factor of two, both for the 
differential distribution and the inclusive Z/"f* cross sec- 
tion used for normalization. A sample of events is then 
generated with ALPGEN v2.13, again with up to three par- 
tons in the matrix element calculation. The factorization 
scale is set to the sum in quadrature of the mass and pr 
of the Z/"/*, and the renormalization scale set according 
to the CKKW prescription. Parton jets from the matrix 
element calculation are required to have pt> 13 GeV, 
and Ai?(jet, jet) > 0.4. These events are hadronized 
in three ways. First, using herwig v6.510 (with jimmy 
v4.31 35] for multiple parton interactions). Then using 
pythia v6.421 with underlying event tune QW (using 
the <5 2 -ordered shower) and the 2-loop prescription for 
a s . Finally, using pythia v6.421 with the p-r-ordered 
shower [361 ] . for which there is currently no tune usinp 
the CTEQ6.1M PDFs, so instead the Perugia* tune 
using the MRST2007 modified LO (LO*) PDFs 
used. This results in three different ALPGEN predictions, 
and in each case the default matching procedure is ap- 
plied after hadronization, with each parton jet required 
to match a particle level jet with pt > 18 GeV, by re- 
quiring Ai?(jet, jet) < 0.4. To determine the impact of 
the matching to the ALPGEN matrix elements calculation, 
herwig and pythia are also run stand-alone to produce 
three inclusive Z/j* — ► fi + \i~ samples. First using her- 
wig v6.510 with jimmy v4.31 for multiple parton inter- 
actions. Then two using pythia v6.421: one with tune 
QW and the 2-loop prescription for a s ; the other with 
the Perugia* tune and the MRST2007 LO* PDFs. In 
both PYTHIA and herwig, the renormalization and fac- 
torization scales for the hard scatter are set to the mass 
of the Z/'y*, and for the initial and final state showers are 
determined dynamically. For all generators, the particle 
level quantities are extracted as defined earlier. 

The normalized differential cross sections are available 
in Ref. [3§|, and presented binned in A</>(Z,jet) (Figs. 
□ and! Tables U and P), |Ay(Z,jet)| (Figs. 1 and H 
Tables ED and E) and \y bo0 Bt(Z + jet)j (Figs. [5] and El 
Tables fVl and IVI j) . The data points are placed at the bin 
average, defined as point where the differential cross sec- 
tion within the bin, taken from simulation re-weighted 
to match the shape in data, is equal to the measured 
value in the bin [40j • The data are shown with statistical 
uncertainties (inner error bar) and combined statistical 
and systematic uncertainties (outer error bar). For clar- 
ity, only the predictions of NLO pQCD and SHERPA are 
shown with the data in part (a) of each figure. In the 
other parts of each figure, ratios are shown, where the 
data and all other theory predictions are divided through 
by the prediction from SHERPA. We choose to show 
the data only in one ratio, to avoid repeating the data 
uncertainties and statistical fluctuations several times. 



SHERPA is chosen as the common denominator for all ra- 
tios as it provides the best description of the shape of the 
data in most distributions, simplifying the determination 
of trends in other theoretical predictions relative to the 
data. The SHERPA scale uncertainty is shown as a shaded 
band around unity. 

The integrated cross sections are also extracted: 
CTz+jct/cz for the stated Z/'y* and jet selections. These 
are measured to be [122 ± 2(stat.) ± 4(syst.)] x 10~ 3 
for pf > 25 GeV, and [47 ± l(stat) ± 2(syst)] x 10~ 3 
for pf > 45 GeV. In both cases, the pf requirement 
is only made on the Z/7*+jet selection, not the in- 
clusive Z/"f* selection. The corresponding results from 
pQCD are [111 ± 6(scale) ± 2(PDF)] x 10" 3 at NLO 
and [112 ± 20(scale) ± l(PDF)] x 10~ 3 at LO for pf > 
25 GeV, and [40 ± 3(scale) ± l(PDF)] x 10~ 3 at NLO and 
[40 ± 8(scale) ± l(PDF)] x 10~ 3 at LO for pf > 45 GeV. 



TABLE I: The measured cross section in bins of A</>(Z, jet) 
for Z/^y* +jet+V events with pf > 25 GeV, normalized to the 
measured Z/'y* cross section. 





(A0 


1/ux da/dA<j> 


<5cr s tat. 


ScTtotal 






(rad) 


(rad) 


(1/rad) 


(%) 


(%) 






0.0-1.5 


1.09 


2.82 x 10" 4 


12 


+24 


-26 




1.5-2.2 


1.95 


4.22 x 10~ 3 


6.9 


+10 


-11 




2.2-2.5 


2.38 


1.93 x 10" 2 


5.5 


+7.3 


-7 


9 


2.5-2.7 


2.61 


5.27 x 10~ 2 


4.1 


+6.1 


-6 


2 


2.7-2.9 


2.81 


1.13 x 10" 1 


2.8 


+4.5 


-4 


6 


2.9-tt 


3.04 


3.32 x 10" 1 


1.7 


+3.6 


-3 


3 



TABLE II: The measured cross section in bins of A(j>(Z, jet) 
for Z/'y* +jet+V events with pf > 45 GeV, normalized to the 
measured Z/'y* cross section. 



Acj> 


(A<^> 


l/a x dcr/dA4> 


<5cr s tat. 


<5<Ttotal 






(rad) 


(rad) 


(1/rad) 


(%) 


(%) 






0.0-1.5 


1.09 


4.88 x 10" a 


26 


+42 


-44 




1.5-2.2 


1.95 


7.61 x 10" 4 


16 


+18 


-18 




2.2-2.5 


2.38 


6.57 x 10~ 3 


9.1 


+11 


-11 




2.5-2.7 


2.61 


1.91 x 10" 2 


6.9 


+8.2 


-8 


6 


2.7-2.9 


2.81 


3.83 x 10~ 2 


4.8 


+6.0 


-6 


5 


2.9-tt 


3.04 


1.35 x 10 _1 


2.5 


+3.5 


-3 


5 



Where it is valid, the NLO pQCD calculation provides 
a good description of the data and is a significant im- 
provement in both shape and uncertainty over LO. How- 
ever, an overall normalization difference of just over 1 
standard deviation of the combined data and NLO the- 
oretical uncertainties is observed, and slightly larger in 
the pf > 45 GeV sample. Of the event generators, her- 
wig shows significant disagreement with data in both 
A<p(Z, jet) and \Ay(Z,jet)\. The trend in the |Ay(Z,jet)| 
is consistent with the description of the leading jet rapid- 
ity Q by herwig, and is significantly improved when in- 
terfaced to the ALPGEN matrix element calculation. The 
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FIG. 1: The measured normalized cross section in bins of A(f>(Z, jet) for Z/~f* +jet+X events for pf> 25 GeV. The distribution 
is shown in (a) and compared to fixed order calculations in (b), parton shower generators in (c), and the same parton shower 
generators matched to ALPGEN matrix elements in (d). All ratios in (b), (c), and (d) are shown relative to SHERPA, which 
provides the best description of data overall. 



TABLE III: The measured cross section in bins of \Ay(Z,]et)\ TABLE IV: The measured cross section in bins of \Ay(Z, jet) | 
for Z/y*+jet+X events with Pt> 25 GeV, normalized to the for Z/j* +jet+X events with Pt> 45 GeV, normalized to the 
measured Z/"/* cross section. measured Z/'y* cross section. 
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predictions of the new Perugia* tune of pythia pro- 
vide a good description of the data in |Ay(Z,jet)| and 
|yboost(Z + jet) I , but not in A0(Z,jet). The modelling 
of A</)(Z,jet) is improved when pythia is interfaced 
to ALPGEN. In general, the three predictions obtained 
from ALPGEN provide a good description of the shape 
of \Ay(Z,iet)\ and |y b oost(^ + jet)| for p|> 45 GeV, 
but perform less well for the other distributions mea- 



sured. Overall, SHERPA provides the best description of 
the shape of data, but shows a significant normalization 
difference. Further, in the sample with > 45 GeV the 
SHERPA description of \yboost(Z + jet) | shows a slope rela- 
tive to the data, which may also be present in \Ay(Z, jet) | 
though it is less clear. All event generators suffer from 
significant scale uncertainties, of comparable size to the 
uncertainty on the LO pQCD prediction. For | Ay(Z, jet) | 
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FIG. 2: The measured normalized cross section in bins of A(f>(Z, jet) for Z/~f* +jet+X events for pf> 45 GeV. The distribution 
is shown in (a) and compared to fixed order calculations in (b), parton shower generators in (c), and the same parton shower 
generators matched to ALPGEN matrix elements in (d). All ratios in (b), (c), and (d) are shown relative to SHERPA, which 
provides the best description of data overall. 



TABLE V: The measured cross section in bins of 
\yboost(Z + jet) | (denoted \yh\) for Z/7*+jet+X events with 
Px> 25 GeV, normalized to the measured Z/-y* cross section. 



TABLE VI: The measured cross section in bins of 
|2/boost(Z + jet) | (denoted |i/b|) for Z/7*+jet+X events with 
Pt> 45 GeV, normalized to the measured Z/-y* cross section. 
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and \yboost{Z + jet)|, these uncertainties change the nor- 
malization with little effect on the shape. In the case 
of A(/>(Z,jet), there is some shape dependence, indicat- 
ing that ALPGEN could be brought into better agreement 
with the data through a lower scale choice. Studies car- 
ried out by shifting the renormalization and factorization 
scales down in ALPGEN and the corresponding pythia 



showering confirm that this is the case. 

In summary, we have presented the first measurements 
at a hadron collider of the Z/j*+jct+X normalized dif- 
ferential cross section in A^(Z,jet), \Ay(Z, jet) | , and 
|yboost(-Z + jet) | . The measurements were made using a 
sample corresponding to 0.97 ± 0.06 fb -1 of integrated 
luminosity recorded by the DO experiment in pp colli- 
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sions at ^/s = 1.96 TeV. These measurement test the 
current best predictions for vector boson + jet produc- 
tion at hadron colliders, and are essential inputs for the 
tuning of event generators. Improving the modeling of 
this important signal will lead to increased sensitivity of 
searches for rare and new physics. 
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TABLES OF RESULTS 



Tables I VIII to IXIII contain the normalized differen- 
tial cross sections, l/a(Z/j*) x da(Z/j* +jet)/dX, for 
the three angular variables A0(Z,jet), |Ay(Z,jet)| and 
|2/boost(Z + jet) | . For each bin we present the bin edges, 
bin average (defined in the main text), measured value, 
statistical uncertainty, uncorrelated systematic uncer- 
tainty, and the systematic uncertainties from the follow- 
ing sources, which are correlated across all bins and all 
distributions: 

1. A0, Ay, j/boost re-weighting; 

2. Jet y re-weighting; 



3. Z/"f* px re-weighting; 

4. Jet px re-weighting; 

5. Data overlay on simulation; 

6. Muon 1/pt resolution; 

7. Muon identification efficiency; 

8. Jet identification efficiency; 

9. Jet energy resolution; 
10. Jet energy scale. 
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TABLE VII: The measured cross section in bins of A<j>(Z,jet) for Z/-y* +jet+X events with pf > 25 GeV, normalized to the 
measured Z/7* cross section. 
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TABLE VIII: The measured cross section in bins of \Ay(Z, jot) | for Z/y*+jet+X events with pf> 25 GeV, normalized to the 
measured Z/j* cross section. 
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TABLE IX: The measured cross section in bins of \yboost(Z + jet) (\yb\) for Z/j'+jet+X events with pf > 25 GeV, normalized 
to the measured Z/7* cross section. 
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TABLE X: The measured cross section in bins of Ac/>(i?, jet) for Z/y*+jet+X events with pf > 45 GeV, normalized to the 
measured Z/7* cross section. 
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TABLE XI: The measured cross section in bins of |Aj/(Z,jet)| for Z/^+jet+X events with pf > 45 GeV, normalized to the 
measured Z/7* cross section. 
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TABLE XII: The measured cross section in bins of \yboost(Z + jet) | (|j/b|) for Z/y* +jet+X events with pf > 45 GeV, normalized 
to the measured Z/7* cross section. 
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